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We calculate, within the single-loop or equivalently the Hartree-Fock approximation (HFA), the 
finite-temperature interacting compressibility for three-dimensional (3D) Dirac materials and renor¬ 
malized quasiparticle velocities for 3D and two-dimensional (2D) Dirac materials. We find that in 
the extrinsic (i.e., doped) system, the inverse compressibility (incompressibility) and renormalized 
quasiparticle velocity at fc = 0 show nonmonotonic dependences on temperature. At low temper¬ 
atures the incompressibility initially decreases to a shallow minimum with a In T dependence. 
As the temperature increases further, the incompressibility rises to a maximum and beyond that it 
decreases with increasing temperature. On the other hand, the renormalized quasiparticle velocity 
at fc = 0 for both 2D and 3D Dirac materials first increases with T^, rises to a maximum, and 
after reaching the maximum it decreases with increasing temperature. We also find that within the 
HFA, the leading-order temperature correction to the low-temperature renormalized extrinsic Fermi 
velocity for both 2D and 3D doped Dirac materials is ln(l/r). 

PACS numbers: 71.55.Ak, 71.18.-|-y 


I. INTRODUCTION 

Recent experimental discovery of three-dimensional 
(3D) Dirac semimetalJi®^ has opened up a new avenue 
for realizing exotic quantum phenomena such as qu antum 
anomalous Hall effecl!i^, giant diama gnetisrrP^R^, oscil¬ 
lating qu antum spin Hall effecP^I^Il^ nonlocal quantum 
oscillatiorP^Ull and negative or linear quantum magne- 
toresistanc^2lHlll This novel material is characterized by 
linear energy band dispersion in the 3D momenta near 
the Dirac point (i.e., the touching point between valence 
and conduction bands). This Dirac cone structure has 
been experimentally observed using ARPES measure¬ 
ment in sev eral materials such asNasBfiH^, Cd 3 As 2 p^, 
TaA^SMUJi^ NbAji^, and Taff^^HH jf the degeneracy 
of the Dirac cones is lifted by breaking either the time- 
reversal or inversion symme try, t hen there are nontriv¬ 
ial Fermi arc surface state^^MMI] unusual transport 
propertie^^ in the resultant Weyl materials. (Since the 
theory developed in this paper applies to Dirac and Weyl 
materials with the appropriate numerical degeneracy fac¬ 
tor correctly incorporated, in this current work we do 
not distinguish between them and hereafter we will refer 
to both as 3D Dirac materials, systems or semimetals). 
In addition to transport properties, the thermodynamic 
property such as electronic compressibility is also of great 
interest as it provides information about the ground state 
many-body effects in interacting Dirac systems with chi¬ 
ral linear band dispersion. The interacting electronic 
compressibility has been extensively studied both theo¬ 
retically and experimentally inparabolic-E^^^ and linear- 
dispersion band materialpSHSU particular, the inter¬ 
acting compressibility in 3D metal^^, two-dimensional 
(2D) semiconductor structure^^HSIl^ a^d graphen^^sEH 
has often been the central quantity in comparing theory 
and experiment in the study of many-body ground state 


properties of strongly interacting electron systems. 

In this paper, we develop a theory for the finite- 
temperature compressibility of a 3D Dirac material and 
renormalized quasiparticle velocities of 2D and 3D Dirac 
materials within the single-loop Hartree-Fock approxima¬ 
tion (HFA). We also give the theoretical results for the 
zero-temperature Hartree-Fock (HF) and random phase 
approximation (RPA) compressibilities in Appendices, 
comparing the two approximations. We note that al¬ 
though we refer to our theory as the HFA (the preferred 
terminology in condensed matter physics), the theory 
is the so-called single-loop field theoretic approximation 
where the calculation is carried out to the leading order 
in the coupling constant in a perturbative sense. The 
leading-order nature of the single-loop HFA also implies 
that the temperature-dependence we obtain would be 
correct to the leading order in the interaction coupling 
constant. We emphasize that the HFA is the leading- 
order perturbative expansion in the bare Coulomb in¬ 
teraction, which is known to produce the exact result for 
the interacting compressibility to the leading order in the 
dimensionless coupling constant (i.e., the effective fine 
structure constant for the Dirac material). We further 
remark that for small coupling constant, the compress¬ 
ibilities calculated within the HFA and RPA have exactly 
the same leading-order divergent term as discussed in the 
Appendix where we present our RPA results in terms of 
an expansion in the dynamically screened Coulomb in¬ 
teraction. 

For a 3D Dirac semimetal, the effective noninteracting 
low-energy Hamiltonian around a Weyl node (which for 
our purpose, leaving out any real spin degeneracy, is also 
a Dirac point where the linear conduction and valence 
bands meet) is 


Ho = • cr, 


( 1 ) 
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where DiJ’ is the noninteracting Fermi velocity (which is 
assumed to be momentum independent), k is the momen¬ 
tum and (T is vector of the Pauli operator. Throughout 
this paper, we set S = 1. The eigenenergies of the Hamil¬ 
tonian are given by es{k) = svpk with the eigenstates 
(for k, ^ op 



cosds/2 \ 
*‘^sind^/2j ’ 


( 2 ) 


where s = ± denotes the conduction and valence bands, 
respectively, ip = aictan{ky/kx), cos'd+ = A: 2 /|k| and 
= TT — i?_|_. For undoped or intrinsic Dirac semimetal, 
the chemical potential is zero (at the Weyl node) while 
the doped or extrinsic Dirac semimetal has finite chemical 
potential. The extrinsic case is generic as the chemical 
potential can be easily shifted away from the Weyl node 
by disorder (and even for a nominally undoped system, 
the chemical potential is located typically either in the 
conduction or the valence band). Moreover, the chemical 
potential can also be adjusted experimentally via surface 
doping. Throughout this paper, we focus on extrinsic 
Dirac materials with positive Fermi energy {Ep >0). At 
zero temperature, the noninteracting chemical potential 
of a 3D Dirac semimetal is the same as the Fermi energy, 
i.e.. 


Mo = Ep = vpkp, 


where the Fermi wave vector is given by 


kp = 



( 3 ) 

( 4 ) 


with n and g being the free carrier density and degener¬ 
acy factor (number of distinct valleys, spin, etc.), respec¬ 
tively. 

The compressibility is defined as the change of vol¬ 
ume with pressure K = —V~^{dV/dP)T,N where K 
is the compressibility, V is the volume, P is the pres¬ 
sure, N is the particle number, and T is the te mpera- 
ture of the system. It can be shown thalP^^^ the in¬ 
verse compressibility (incompressibility) is related to the 
change of chemical potential with respect to carrier den¬ 
sity by K~^ = n^dp/dn. In 2D electron gas, the quantity 
dfx/dn can be obtain ed expe rimentally by measuring the 
quantum capacitanc^^^^^SElHlll or directly measuring the 
density-dependent che mical p otential using a scanning 
single-electron transistoM^^^. At zero temperature, the 
noninteracting incompressibility is simply 


1 2 

Kq dn Dq ’ 


( 5 ) 


/g 2\1/3 

where Dq = ;^ ( j is the density of states at the 
Fermi surface of a noninteracting 3D Dirac semimetal. 
For a 3D Dirac semimetal, dpo/dn oc due to the 

linear energy dispersion while for a parabolic-band 3D 
electron gas (3DEG), dpo/dn oc 


Interaction among the electrons profoundly affects the 
compressibility since the chemical potential is renormal¬ 
ized by many-body interaction effects. In particular, it 
is known that interaction could drive the electronic com¬ 
pressibility negative in 3D and 2D parabo lic-ba nd elec¬ 
tron gases (as elucidated in Ref.^^ and Refs.l^SlMl, respec¬ 
tively) whereas in graphene, the prototypical 2D Dirac 
material, electron-electron interaction renormalizes the 
chemical potential, but does not lead to a negative com- 
pressibilitjf^. In the current work, we address the subject 
of electron-electron interaction effects on the compress¬ 
ibility of 3D Dirac materials. In the usual 3D metals, 
the electronic temperature scale is Tp ^ 10^ K (where 
Tp = Ep/kp is the Fermi temperature with kp being the 
Boltzmann constant), and hence the temperature depen¬ 
dence of many-body renormalization effects is not a par¬ 
ticularly relevant topic of research in 3D metals although 
for lightly doped 3D semiconductors such temperature 
dependence may become relevant. In 2D systems, how¬ 
ever, both for 2D semiconductor structures and graphene, 
the chemical potential (i.e., the 2D carrier density) can be 
tuned almost arbitrarily using an external gate voltage, 
and therefore, temperature-dependent many-body effects 
take on considerable significance. In the current work, we 
consider interaction effects for both the zero-temperature 
and finite-temperature compressibilities in 3D Dirac ma¬ 
terials since the chemical potential can be tuned in 3D 
Dirac systems by extrinsic doping and for low doping, 
one expects the temperature dependence of the many- 
body renormalization effect to be significant by virtue of 
the Fermi energy being small. In this work, we assume 
that the bare electron-electron interaction in the system 
is the long-range Coulomb interaction as appropriate for 
electrons in standard semiconductors, semimetals, and 
metals. 


The paper is structured as follows. In Sec. [ITj we dis¬ 
cuss the finite-temperature compressibility of 3D Dirac 
semimetals within the HFA. More specifically, we give 
the full numerical results of the finite-temperature com¬ 
pressibility and its analytic expression at low temper¬ 
atures. In Sec. |III[ we give the full numerical results 
and low-temperature analytic expression of the finite- 
temperature HF renormalized quasiparticle velocities at 
k = 0 and k = kp for 2D and 3D Dirac materials. A 
summary is given in Sec. IV In the appendices, we give 
detailed calculations of the zero-temperature HF com¬ 
pressibility of 3D Dirac materials (Appendix [A|) as well 
as the zero-temperature RPA compressibility of 2D and 
3D Dirac materials (Appendix [ b|). We note that the zero- 
temperature theoretical results are provided in the Ap¬ 
pendices with the main text focusing on the theoretically 
challenging finite-temperature results. We believe that 
providing both T = 0 and finite-T results in the same 
paper elucidate different aspects (e.g., temperature and 
density dependence) of the many-body renormalization 
effects in Dirac materials within one coherent formalism. 
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II. FINITE-TEMPERATURE 
COMPRESSIBILITY 


In this section, we present the finite-temperature com¬ 
pressibility of 3D Dirac materials in the single-loop (or 
the leading-order in interaction) theory, which is equiv¬ 
alent to taking into account the effect of exchange- 
interaction. More specifically, we focus on the extrin¬ 
sic or doped Dirac semimetal with Ep > 0. The zero- 
temperature HF compressibility is given in Appendix [A| 
The noninteracting extrinsic chemical potential at fi¬ 
nite temperature is determined by the conservation of 
total electron density: 


V 


^4+)(k) + 


V 


[n'p ^(k) 



= n, 


( 6 ) 


where n'p'^ is the Fermi-Dirac distribution at temper¬ 
ature T with energy £± = ±vpk and n = gkp/{6Tr'^) 
is the zero-temperature carrier density due to doping. 
From Eq. (§, we obtain the self-consistent equation for 
the finite-temperature noninteracting chemical potential 
Po as 


/ip 

Ep 


{ntr 



( 7 ) 


where t = T/Tp is the normalized temperature with 
Tp = Ep/kp being the Fermi temperature. Solving 
Eq. 0, we obtain the noninteracting chemical potential 
in the limit of low and high temperatures a^^ 


Within the HFA, the interacting chemical potential is 
the sum of kinetic energy and exchange self-energy. The 
exchange self-energy is given by: 

^Tik) =-T.j ^n^(6'q)Ek-qE,,,(k,q), (10) 

where s,s' = ± are the band indices, Uq = Ape^Khcq^) 
is the “bare” long-range (^ 1/r) Coulomb interaction in 
the momentum space with k being the background di¬ 
electric constant of the material, ^s'q = s'vpq — go is the 
free dispersion, np{^s'ci) is the Fermi-Dirac distribution, 
andFss/(k, q) = |(ks|qs')| = (l-|-ss'cosd )/2 is the over¬ 
lap of the eigenstates with 6 being the angle between k 
and q. We note that the Hartree part of the self-energy 
vanishes exactly for the homogeneous system since it is 
canceled by the positive background charge necessary for 
neutrality, thus leaving the exchange or the Fock self¬ 
energy as the only leading-order single-loop term. The 
exchange self-energy can be split into the intrinsic self¬ 
energy (S™*), which arises from the contribution at zero 
chemical potential and the extrinsic self-energy (S®^*), 
which takes into account the correction to the self-energy 
due to nonzero chemical potential: 

K\k) = -J ^Vk-qE,_(k,q), ( 11 a) 

Sr(fc) = -E / ^^n^(e.'q)Ek-qE«,,(k,q), 

(lib) 


Ep ] I ..... 


The incompressibility is given by 1/K = n^dg/dn. Since 
I/Kq = n^dEp/dn, we have Kq/K = dg/dEp = 
^ finite-temperature noninteracting incom¬ 

pressibility is then given by the following analytical for¬ 
mulas for r <C Tp and T '^Tp, respectively. 


where 5np[^s'<^ = np{^s'q) — (1 — s')/2 is the differ¬ 
ence between the electron occupation number of finite- 
temperature extrinsic 3D Dirac semimetals and the in¬ 
trinsic T = 0 case. 

Since the chemical potential (consequently, the com¬ 
pressibility) is only affected by the electron energy at 
Fermi momentum kp, in the following, we will evaluate 
the self-energy at k = kp. Considering extrinsic 3D Dirac 
materials with Ep > 0 (s = -1-1), we have the intrinsic 
self-energy a^^ 


Kq 

KUT) 


3 ’ 

for t ^ 1 , 

3 

772 ^ 2 ’ 

for t 1 , 


( 9 ) 


where Kq and Kni(T) are the zero-temperature and 
finite-temperature noninteracting compressibility, re¬ 
spectively. Thus, to the leading-order in t, the in¬ 
compressibility Kni(T)~^ increases quadratically and de¬ 
creases inverse quadratically with temperature in the 
limit of low and high temperatures, respectively. This 
is to be contrasted with the behavior of noninteracting 
incompressibility of a parabolic 3DEG (i.e., normal 3D 
metals and doped semiconductors), where it increases 
quadratically and linearly with temperature in the limit 
of low and high temperatures, respectivelji^. 


^ikp)=- 


avpkp f 1 
2r] 

3(ry — rf — rf) — 1 


6ry2 3 


1-4 

rjz 


+ - 


12773 


In 


1 -I- 77 


1 — 77 


( 12 ) 


where a = e^KKVp) and 77 = kp/kc with fc. being the 
ultraviolet momentum scale beyond which the linear dis¬ 
persion no longer holds. We emphasize that fc. is the 
ultraviolet momentum scale necessary for cutting off the 
high-momentum logarithmic divergence associated with 
the linear Dirac spectrum. Eor 3D Dirac systems, fc. 
should typically be of the order of the inverse of the lat¬ 
tice length scale (or perhaps even a much smaller mo¬ 
mentum scale in the realistic material where the band 
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energy dispersion starts deviating from linearity). Us¬ 
ing the renormalization group (RG) theory it is possible 
to express the interacting compressibility at one density 
(i.e., one kp value) simply in terms of that at another 
density, thus completely eliminating the unknown ultra¬ 
violet scale associated with kc as has been extensively 
discussed recently in Ref.^^. Note that the interaction 
coupling constant a = /{kvf) is the ratio of the in¬ 
teraction to kinetic energy (the effective hne-structure 
constant), which is density independent for Dirac mate¬ 
rials. For the case of parabolic dispersion, the coupling 
is density dependent and typically denoted by Vg where 
Tg oc for regular 3D parabolic-band systems. The 

extrinsic self-energy for low temperatures {t ^ 1 ) can be 
Sommerfeld expanded, yielding 

ET'*(A:f) = - (2 -I- 41n2 - Tr^t^lnt) , 

(13) 

where A = 2^ + ln(32/7r) — 121nG ~ 0.49 with 7 ~ 
0.577 being the Euler’s constant and G ~ 1.282 being 
the Glaisher’s constant. 

Within the HFA, the interacting chemical potential is 
the sum of noninteracting part and exchange self-energy 
at Fermi momentum: 


/i = /ro + S“(fcF) + S;"‘(A:F). (14) 


Differentiating the chemical potential with respect to the 
density, we then obtain the low-temperature interacting 
incompressibility Kq/K{T) as 


= 1 

K{T) 3 


a 

TT 


9 

2+2 


kp 

kr 


1 2 , „ 

-I- 1 - - In 2 

3 3 


6 


(15) 


where 


giji) = 


127?3 


65 ^ — 4:7] — Aff In 

1-4 

- 



+(2 — 3t] — 3 ? 7 ^) In 


l + rj 


1-77 


= ^ + ^ln?? + G(?7)- 


(16) 


Similar to the case of graphene and parabolic-band sys- 
tem^^, the low-temperature HF incompressibility of a 
3D Dirac material shows a leading-order In t temper¬ 

ature dependence followed by a next-to-leading-order 
term. 

Figure [ 1 ] shows the plots of the exact numerically cal¬ 
culated single-loop finite-temperature incompressibility 
Kq/K. The results are computed by using wave vector 
cutoff kc = 27r/ac (oc = 16 A is the geometric mean of 
the lattice constants in Gd 3 As 2 ), degeneracy 5 = 4 and 
Fermi velocity vp = 10® cm/s. We provide our numerical 
results using the parameters for Gd 3 As 2 unless otherwise 






FIG. 1. (Color online) (Top) Incompressibility KojK as a 
function of normalized temperature t = T jTp for three differ¬ 
ent carrier densities n (expressed in the unit of 10^® cm“®) and 
dielectric constants: (a) k = 1 and (b) 36. Inset of (a) shows 
the low-temperature behavior of the incompressibility. Solid 
lines denote the numerical results and red dashed lines denote 
the asymptotic low-temperature analytical results. (Middle) 
Incompressibility Kq/K as a function of carrier density n for 
different values of normalized temperature t and dielectric 
constants: (c) «; = 1 and (d) 36. (Bottom) Semilog plot of 
Kq jK as a function of carrier density n for different values of 
temperature T and dielectric constant: (e) k = 1 and (f) 36. 
The above plots are calculated using the parameters: wave 
vector cutoff {Uc = 16 A), degeneracy 5 = 4, and 

Fermi velocity uf = 10® cm/s. 


stated. Figs, [^a) and (b) show the temperature depen¬ 
dence of the incompressibility Kq/K for three different 
values of carrier density n = 0 . 1 , 1 and 10 (given in the 
unit of 10 ^® cm“®) with dielectric constants: k = 1 and 
K = 36 (dielectric constant of Cd 3 As 2 at 4.2 K)P. As can 
be seen from the hgures, the temperature dependence is 
nonmonotonic. At “low” temperatures (i.e., t <C 1) cor¬ 
responding to high values of doping density (and conse- 
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quently Tp can be quite large), the leading-order tem¬ 
perature dependence is In t where it develops a shallow 
minimum with the size of the minimum proportional to 
the value of a. Due to the smallness of the prefactor air/6 
[Eq. @], this minimum is very shallow as shown in the 
inset of Fig. [^a). The low-temperature asymptotic an¬ 
alytical results [Eq. (15l], shown as red dashed lines in 
Fig-Qa), agree precisely with the numerical results. As 
the temperature increases further, the incompressibility 
begins to increase until it reaches a maximum and beyond 
that it decreases with increasing temperature. Unlike the 
parabolic-band electron gas where the exchange-energy 
contribution to the compressibility vanishes at high tem¬ 
perature^^, for Dirac materials there is a finite contri¬ 
bution of the exchange self-energy to the compressibil¬ 
ity at large t. This interesting persistence of quantum 
many-body effects into the classical T ^ Tp regime pre¬ 
sumably arises from the fact, pointed out in Ref.l^, that 
Dirac systems with linear energy dispersion are mani¬ 
festly nonclassical and do not have any classical analogs 
since a linear-in-momentum energy dispersion is mani¬ 
festly quantum-mechanical and not allowed classically. 
We note that a large dielectric constant suppresses the 
electron-electron interaction strength (since the coupling 
constant is inversely proportional to the background di¬ 
electric constant) which weakens the density dependence 
of Kq/K as can be seen by comparing Figs.[^a) and (b). 



FIG. 2. (Color online) (a) Incompressibility KojK as a func¬ 
tion of normalized temperature t = T/Tp for three different 
values of dielectric constant k = 1, 3 and 36. The black solid 
line (NI) denotes the noninteracting incompressibility, (b) 
Incompressibility Kq /A as a function of dielectric constant n 
for different values of normalized temperature t. The plots 
are calculated using the parameters: n = 10^® cm“®, gi = 4, 
Vp ~ 10® cm/s and A = 2p/a^ (a^ = 16 A). 

The middle panel of Fig. shows the density depen¬ 
dence of Kq/K for several fixed values of t. For fixed t, 
the incompressibility decreases with the carrier density n 
due to the decrease in the ratio of kc to kp (where kc is 
assumed to be fixed for a given material and kp increases 
with density n). This is to be contrasted with the case 
of regular 3D parabolic-band system^P^ where for fixed t, 
the incompressibility increases with n. Note that a large 
background dielectric constant suppresses the density de- 



FIG. 3. (Color online) (a) Incompressibility KojK as a func¬ 
tion of wave vector cutoff kc/kp for different temperatures 
t = T/Tp. (b) Incompressibility KojK as a function of nor¬ 
malized temperature t for different values of wave vector cutoff 
kc/kp. The above plots are calculated using a = 2.2. 


pendence of the incompressibility. The nonmonotonic 
dependence of the incompressibility on temperature can 
also be seen from Figs, [^c) and (d). The bottom panel 
shows the semilog plot of the incompressibility Kq /K as 
a function of carrier density n. For a fixed temperature 
T, at a high carrier density which corresponds to the 
low-t limit, Kq/K is a decreasing function of n mani¬ 
festing density dependence while at a low carrier 

density corresponding to the high-t limit, Kq/K increases 
as ri^/^. This dependence on n (where n (x t~^) follows 
mainly from the temperature dependence of the nonin¬ 
teracting part [Eq. 

Figure [^a) displays the temperature dependence of 
Kq/K for different values of the background dielectric 
constant k. The noninteracting part of the compressibil¬ 
ity becomes more dominant for larger dielectric constants 
as the electron-electron interaction strength is suppressed 
in the presence of large dielectric constant. The ultra¬ 
violet momentum cutoff dependence of Kq/K is shown 
in Fig. As shown in the figure, the incompressibil¬ 
ity increases with increasing momentum cutoff [with a 
\n{kc/kp) dependence for large cutoff]. To better aid ex¬ 
perimentalists in interpreting the cutoff dependence, in 
Fig. we plot Ao(n)/A(n) for different values of densi¬ 
ties n against KQ{nQ)/K( uq) at no = 10^® cm“® as the 
momentum cutoff kc/kp is varied from 1 to 100. 

We note that experimentally the ultraviolet cutoff kc is 
of course unknown (although using a cutoff correspond¬ 
ing to the band momentum where the relevant band dis¬ 
persion deviates from linearity or even the lattice cutoff 
should both be fairly decent approximations), and there¬ 
fore, the interaction coupling should be considered as a 
running coupling in the standard RG sense as discussed 
in details recently in Ref.l^. A practical way to compare 
theory with experiment would be to eliminate kc and 
express the experimental compressibility at one density 
with that at another density. Taking one of these den¬ 
sities to be a fixed fiducial density, the measured com¬ 
pressibility at an arbitrary density could be expressed in 
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Ko(no)IK(no) Ko(no)IK(no) Ko{no)IK{no) 


FIG. 4. (Color online) (Top) Incompressibility Ko{n)/K{n) 
vs Ko(no)/K(no) as ka/kp is varied from 1 to 100 for different 
temperatures T and various values of density n (given in the 
unit of 10^® cm“^): (a) n = 0.01, (b) 0.1, and (c) 10. (Bot¬ 
tom) Incompressibility Koin)/K(n) versus Ko{no) / K{no) as 
kc/kp is varied from 1 to 100 for different densities n (xlO^® 
cm“®) and various values of temperature: (d) T = 300 K, 
(e) 1000 K, and (f) 3000 K. In all of the above plots we set 
no = 10^® cm“^ and a = 2.2. 


terms of this fiducial compressibility and thus compared 
with the theory. Such a comparison, as described already 
in the literatur^SoMl^ directly tests the running coupling 
constant idea of quantum field theoretic RG flow in the 
specific context of solid-state Dirac materials. 


III. FINITE TEMPERATURE RENORMALIZED 
QUASIPARTICLE VELOCITY 

In this section we calculate the finite-temperature 
renormalized quasiparticle velocity for extrinsic 3D and 
2D Dirac materials with Ep > 0 within the single-loop 
HFA. We note that the interaction theory for the com¬ 
pressibility developed in the last section implicitly incor¬ 
porates many-body velocity renormalization effects, but 
it is useful to provide the explicit results for the quasi¬ 
particle velocity renormalization in Dirac materials as a 
function of temperature since the quasiparticle velocity 
enters as a central quantity determining many proper¬ 
ties of Dirac systems. There is no available theory for 
the temperature dependence of the quasiparticle veloc¬ 


ity renormalization in 3D and 2D Dirac materials. Here, 
we provide temperature dependent results for both cases. 
In particular, we give the full numerical results for the 
one-loop HF renormalized quasiparticle velocity at A: = 0 
together with the low-temperature analytic expressions 
for the one-loop HF renormalized quasiparticle velocity 
at fc = 0 and k = kp. Although the bare velocity vp 
in Eq. Q for the effective linear band energy disper¬ 
sion is by definition momentum independent for Dirac 
systems, interaction-induced many-body effects may in¬ 
troduce nonlinearity or momentum dependence since the 
self-energy is in general momentum dependent. In Ap¬ 
pendix 1^ we present the detailed calculation of zero- 
temperature HF quasiparticle velocity for 3D Dirac ma¬ 
terials. Within the HFA, the renormalized quasiparticle 
velocity for the doped (extrinsic) system with Ep > 0 is 
defined by 


^ext(A^) — 


d[e+{k') + E“)*(fc') -b S^’^t(fc')] 


dk' 


(17) 


k'=k 


A. 3D Dirac material 


At low temperatures, the leading-order temperature 
correction to the renormalized HF quasiparticle velocity 
of a 3D Dirac semimetal at A; = 0 can be calculated to 
be [using Eq. @] 


Vext{k = 0) =Vp <1 + 


a 

Stt 


2 In 


I ^2.2 

+ TT t 


■ (18) 


Thus the low-temperature renormalized quasiparticle ve¬ 
locity at A; = 0 has a leading-order temperature depen¬ 
dence without any logarithmic term in temperature. Fig- 
ure[^a) shows the full numerical results of Uext(A: = 0) for 
3D Dirac materials. As shown in the figure, Uext(fc = 0) 
has a nonmonotonic dependence with temperature where 
it first increases with temperature, reaches a maximum 
and then decreases with increasing temperature t. The 
low-temperature analytical results [Eq. (18l], shown as 
red dashed lines in the figure, are in a perfect agreement 
with the numerical results. 

To the leading-order in A, the renormalized extrinsic 
Fermi velocity [n^ = next(fci^)] at low temperatures is 


< 1 H— 

TT 


tb 


(19) 


where g{r]) is given by Eq. (16) and B = 


ln(3-^/7r^) ~ 1.79. Theln(l/A) dependence of the renor¬ 
malized extrinsic Fermi velocity implies that Vp^ooas 
A —>■ 0. So within the HFA, the renormalized extrin¬ 
sic Fermi velocity diverges logarithmically with temper¬ 
ature. We expect this logarithmic divergence to disap¬ 
pear if the correlation effects are taken into account as, 
for example, in the RPA theory where the perturbative 
expa nsion is carried out in the screened Coulomb inter- 
actiorP^mm This logarithmic zero-temperature diver¬ 
gence of the renormalized Fermi velocity (obtained at 
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/cf) is a well-known pathology of the HFA using the un¬ 
screened Coulomb interaction for the many-body pertur¬ 
bative expansion. This logarithmic HFA-induced Fermi 
surface divergence has b een known for 3D metals for 
more than 50 year^^ZEH^ jl; happens for Dirac ma¬ 
terials only when there is a Fermi surface, i.e., when the 
system is doped. The corresponding T = 0 logarithmic 
divergence in extrinsic graphene has been discussed in 
depth in Ref.l^, and the same description applies here 
for the 3D Dirac system. We emphasize that unlike 
the ultraviolet logarithmic divergence associated with the 
momentum cutoff kc (where the divergence arises from 
kc —t oo) which is a defining property of the linearly 
dispersing Dirac-Weyl systems, the infrared logarithmic 
divergence arises specifically from the long-range nature 
of the electron-electron interaction and is generic to all 
electron liquids with Coulomb interaction. This infrared 
divergence, being a property specific to the Fermi surface, 
does not happen at the band bottom (fc = 0), and there¬ 
fore, there is no such infrared divergence issue for the 
undoped (intrinsic) Dirac system with no Fermi surface. 
Thus, the velocity renormalization at fc = 0 [Eqs. (181] 
does not have the infrared divergence while the renormal¬ 
ized extrinsic Fermi velocity at kp [Eqs. (19)] possesses 
the infrared divergence. We note that, in contrast to the 
compressibility itself, the single-loop velocity renormal¬ 
ization at fc = 0 has a temperature dependence with 
no logarithmic In T leading-order correction. 



FIG. 5. (Color online) Renormalized quasiparticle velocity 
at fc = 0 [Uext(fc = 0)/nF] as a function of rescaled temper¬ 
ature t = T/Tf for (a) 3D Dirac materials with wave vector 
cutoff fcc = 2^1 aa (flc = 16 A) and different carrier densities 
n (xlO^® cm“®) and (b) graphene with wave vector cutoff 
fcc = 27r/ac (oc = 2.46 A) and different carrier densities n 
(xlO^^ cm“^). Solid lines denote numerical results and red 
dashed lines denote the low-temperatnre analytical results. 
The above plots are calculated using a = 2.2 and degeneracy 
ff = 4. 


B. 2D Dirac material 

Although the focus of the current paper is on many- 
body interaction effects in 3D Dirac systems, we pro¬ 


vide in this subsection our FIF single-loop results for 
the temperature-dependent many-body renormalization 
of 2D Dirac quasiparticle velocity, which is the finite- 
temperature generalization of the e xistin g T — 0 velocity 
renormalization results for graphen^^^MI refer to the 
2D Dirac system as graphene with no loss of generality 
because of the vast existing literature on graphend^THIoi 
Solving Eq. (§, we obtain the noninteracting chemical 
potential for graphene in the low- and high-temperature 
limits a^^ 


Mo 

Ep 


2^2 


1 - 
1 


TT t 


4t In 2 ’ 


for t <C 1, 
for t > 1. 


( 20 ) 


Using the above chemical potential, we can then cal¬ 
culate the leading-order temperature correction to the 
low-temperature renormalized quasiparticle velocity for 
graphene at fc = 0 to be 


Uext(fc = 0) = t>F 1 + ^ 


1 f 


-b 




■ ( 21 ) 


Similar to the 3D Dirac case [Eq. (18) above], Uext(fc = 0) 
for graphene also increases as E. The full numeri¬ 
cal results for Uext(^ = 0) in graphene are presented 
in Fig. ib). The renormalized quasiparticle velocity 
Vextik = 0) for graphene has a similar nonmonotonic 
temperature dependence to that for 3D Dirac materials 
where it increases as at low temperatures and decreases 
with temperature at high temperatures, reaching a max¬ 
imum at some intermediate value of temperature. 

To the leading order in temperature t, the low- 
temperature renormalized extrinsic Fermi velocity is 


Vp = Vp 


1 + a 




( 22 ) 


where 


piii) = - 

TT 


1 /(m- 1) 


[E{p)-K,{p)]-- 


r /4\ 

1 

In - 


1. \riJ 

~ 2 


-In 


f'dy^ 

lo 

-)1 

M/J 


K,{y)-E{y)--y^ 


Oiv), 


D = 


2p 


3-b2C-81n2-b21n 


(23a) 
0.266, (23b) 


with C « 0.916 being the Catalan’s constant, Ki{x) and 
E{x) being the complete elliptic integrals of the first and 
second kind. Similar to the 3D case, there is a logarithmic 
dependence of the extrinsic Fermi velocity on tempera¬ 
ture within the HFA which is exactly the same infrared 
divergence associated with the long-range Coulomb in¬ 
teraction discussed before. This infrared logarithmic di¬ 
vergence at kp for the extrinsic system will vanish if the 
correlation effects are taken int o acc ount (e.g., RPA) as 
for the case of zero temperatur^SHM! 



























IV. SUMMARY 


Dirac systems as described in this work. 


In conclusion, we have theoretically obtained the finite- 
temperature interacting compressibility for extrinsic 3D 
Dirac materials and renormalized quasiparticle veloci¬ 
ties for extrinsic 2D and 3D Dirac materials within the 
single-loop HFA using the long-range Coulomb interac¬ 
tion between the electrons. We give the numerical results 
of the incompressibility K/Kq and also its analytic ex¬ 
pression at low temperatures. The HF incompressibility 
K/Kq has a nonmonotonic dependence with temperature 
t = T/Tp. At low temperatures, as the temperature in¬ 
creases the incompressibility first decreases slowly with a 

In t dependence to a minimum. After reaching the min¬ 
imum, as the temperature is raised further it rises to a 
maximum and subsequently decreases with temperature 
at high temperatures. We note that our HFA results are 
the same as the single-loop (in the bare interaction) re¬ 
sults, and our results are therefore exact to the leading 
order in the effective coupling constant. We also present 
the numerical results and low-temperature analytic ex¬ 
pression for the finite-temperature renormalized quasi¬ 
particle velocities at A: = 0 [uext(fc = 0)] of 3D and 2D 
Dirac materials. The renormalized quasiparticle velocity 
at A; = 0 also has a nonmonotonic dependence on temper¬ 
ature, i.e., it first increases with temperature and beyond 
a certain value of temperature, it decreases with increas¬ 
ing temperature. At low temperatures the leading-order 
correction to Vext{k = 0) is P. Our theoretical predic¬ 
tions can be directly verified in the experiments through 
either temperature-dependent measurements at a finite 
doping density (with a fixed Tp = Ep/kp) or through 
density-dependent measurements at a finite temperature 
since the relevant dimensionless variable determining the 
temperature dependence is t = T/Tp, where Tp oc 
and Tp oc in 3D and 2D Dirac materials, respec¬ 
tively. 

Besides the finite-temperature results, we also discuss 
the zero-temperature density-dependent HF and RPA 
compressibilities in the Appendices. As shown in Ap¬ 
pendix]^ the HF incompressibilities for 2D and 3D Dirac 
materials increase logarithmically with the momentum 
cutoff. On the other hand, the RPA incompressibility 
(given in Appendix]^ for 3D Dirac materials has a log- 
log dependence on the cutoff while the RPA incompress¬ 
ibility for 2D Dirac materials increases logarithmically 
with the momentum cutoff but with a smaller prefactor 
compared to its HF counterpart. We note that for weak 
coupling, where the single-loop theory is valid, both HFA 
and RPA give quantitatively identical results (except for 
the infrared logarithmic divergence for the renormalized 
Fermi velocity at k = kp which occurs in HFA, but not 
in RPA), and therefore, a systematic experiment to mea¬ 
sure the compressibility in one of the newly discovered 
3D Dirac materials, most of which have rather low val¬ 
ues of interaction coupling constant (typically, a ^ 0.1), 
as a function of density and temperature should lead to 
direct experimental observation of many-body effects in 
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Appendix A: Zero-temperature HF compressibility 


In this Appendix, we give detailed calculations of the 
zero-temperature interacting compressibility within the 
HFA. The zero-temperature intrinsic and extrinsic ex¬ 
change self-energies are given bjE3 




E^ik) 

where 



avpkp 

TT 




(Ala) 

(Alb) 


1 1 — 

^W =2 + ^‘" 


l-\- X 


I — X 


= i-:^ + Oix% 


, / X 1 2;, 

'‘W = -fe + 3‘" 
= Hx(Ul„x 


1-4 




12x2 


In 


1 -|- X 


1 — X 


Oix^ 


(A2a) 


(A2b) 


Note that at T = 0, the HF self-energy for a parabolic- 
band 3DEG is given by -^^/(^)- 


a. Intrinsic (undoped) Dirac material 


Near the Weyl node {k <C kc), the intrinsic self-energy 
can be expanded as 




avpkc 

TT 




. . 

The first term —avpkd'n introduces a trivial shift in the 
energy while the remaining terms renormalize the quasi¬ 
particle velocity. The renormalized quasiparticle velocity 
for small k {k ^kc) is given by 




= Vp 



(A4) 


where kc = e~^l^kc ~ 0.846A ;e. This renormalization is 
also obtained in Refs IsiMiZllIZl! The logarithmic diver¬ 
gence of Uint(A:) as A: 0 is found only in the intrinsic 
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case as for the extrinsic case this divergence is canceled 
out by a similar divergence in as discussed below. 


b. Extrinsic (doped) Dirac material 


For k kp, the extrinsic self-energy can be written as 




avpkp 

TT 


1 -I- s 


2k 

3kp 




Note that since the Infc term in —)■ 0) cancels 

the same term in E“*(A: —t 0), the renormalized quasi¬ 
particle velocity in the extrinsic case Vextik = 0) = 
vf[1 -\- ln(A:c//cF)] bas no logarithmic divergence for 
/c —>■ 0. This cancellation of logarithmic divergence is 
also found for the case of graphen^^, i.e., 2D Dirac sys¬ 
tems. 

At the Fermi wave vector kp^ the extrinsic self-energy 
is 


^T\kp) 


avpkp 

7T 



(A6) 


So compared to the self-energy of a parabolic-band 
3DEG, the extrinsic self-energy can be written as 


^Tikp) 



^P\kp), 


(A7) 


where Il^^{kp) = —e^kp/{KTT) is the self-energy of a 
parabolic-band 3DEG at Fermi momentum. Figure 
shows the exchange self-energies as functions of wave 
number k/kp for a 3D Dirac material with Ep > 0. As 
for EP*’, there is a logarithmic divergence in the slope of 

E^xt 

as k ^ kp which manifests as the divergence in the 
renormalized extrinsic Fermi velocity. This is the well- 
known infrared HFA divergence associated with the un¬ 
screened Coulomb interaction generic to all electron liq¬ 
uids, which disappears in higher-order approximations, 
e.g., RPA. However, this singularity does not have any 
pathological effect on dpL/dn. On the other hand, E^* 
has a finite slope at k = kp since for extrinsic 3D Dirac 
material with Ep > 0 there is no Fermi surface at k = kp 
in the valence band. 

For doped 3D Dirac materials with Ep > 0, the inter¬ 
acting chemical potential within the HFA is given by 

ti = fio + ^^\kp) + E‘)((\kp). (A8) 


Taking the derivative of the chemical potential with re¬ 
spect to the free carrier density and ignoring terms of 
order kp/kc, we have 


( 


dpL 

dn 


ext 


dfip 

dn 




1 - 


int 


4 2 2kF 

- -k - In —- 

9 3 fee 

qjnt 1 , 2 
TT ^3 3 



(A9) 



FIG. 6. (Color online) Exchange self-energies for a 3D Dirac 
semimetal (with Ep > 0) and for a parabolic-band 3DEG 
as functions of wave number k/kp. Note that EP’’(fc) = 
E.^f{k) + YE^\k). 


where dtipjdn = vf[2p‘^= 

r!)!l*[27r^/(9n^5)]*'^^, and aint = e^/(Kn“*) = 

fint(A:F)]- The exchange interaction changes dp,/dn from 
its noninteracting value djjLp/dn by — f [| + | ln(?^)] and 
from its intrinsic value {dg,/dn)int by — (^-I-§ ln2). So 

for carrier density n = 5 x 10*® cm“®, degeneracy g = 4, 
and kc = 2p ja^ (fle = 16 A), {dpLjdn)e-^t is increased from 
its bare value by « 0.333q! but decreased from its intrin¬ 
sic value by « —0.253ainf Figure [^a) shows explicitly 
the comparison between the noninteracting, intrinsic and 
extrinsic dg/dn. Using dielectric constant k = 36 (for 
Cd 3 As 2 )P^and vp = 10® cm/s, the value of a can be cal¬ 
culated to be « 0.06 and dg/dn is increased by roughly 
2% over its bare value. The effect of dielectric constant 
K on (i9/r/i9n)ext is shown in Fig. Bb). Larger dielec¬ 
tric constant suppresses the electron-electron interaction 
strength which in turn decreases the value of dfi/dn. For 
the parameters: k = 1, vp = 10® cm/s and n = 10*® 
cm“®, {dpt/dn)ext is increased from its bare value by 
about 70%. 

At zero temperature, the incompressibility of an ex¬ 
trinsic 3D Dirac material can be calculated to be 


^ = 1 -- 
K TT 


aiv) + 3 + 


2 

3 


In 2 


(AlO) 


where g{r]) is given in Eq. ( |16[ ) and g = kp/k^- On 
the other hand, within the HEA, the zero-temperature 
incompressibility for a parabolic-band 3DEG is 


Kp ^ X - — 

Kph - ^ 



(All) 


where = (3/47rn)*/®(TOe^K) oc n *7®. So KqIK'^'° in¬ 
creases with the carrier density n. This is contrary to 
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FIG. 7. (Color online) (a) Plot of dfi/dn as a function of 
carrier density n calculated using the parameters: n = 1, 
vf = 10® cm/s, kc = jac. {ac = 16 A), <7 = 4, and a = 2.2. 
The notations “bare”, “int”, and “ext” denote the noninter¬ 
acting, intrinsic and extrinsic quantities. Inset: Plot of renor¬ 
malized intrinsic velocity at kp and extrinsic incompressibility 
as a function of density n. (b) Semilog plot of (9/i/cIn)ext as a 
function of carrier density n for different dielectric constants 
K. The noninteracting duo/dn is shown as a black solid line 
for a comparison. 


the compressibility of a 3D Dirac material, where the 
compressibility decreases with n. 


Appendix B: Zero-temperature RPA compressibility 


In this Appendix, we calculate the leading-order ultra¬ 
violet divergent piece of the zero-temperature RPA com¬ 
pressibility for 2D and 3D extrinsic Dirac materials. The 
RPA self-energy can be written as 


k) — ^ ^ G{iUJn “f k ) 




q) _ 


F«,,(k,k'), (Bl) 


where /3 = l/{kBT), k' = k -f q and e(iDm,q) = 
1 -1- Pqn(fDm,q) is the frequency-dependent dielectric 
constant with the polarizability given by 


q) 


V 


E 

kss' 


npi^sk) - njp(^s'k') ^ 

n \ c c ^^ ^ 

+ qsk — ^s'k' 


(B2) 

The RPA self-energy is obtained by replacing the “bare” 
interaction by dynamically screened interaction (i.e., ge¬ 
ometric ser ies over polarization bubbles with the “bare” 
interaction)P35lIlll, The sum of the series removes the in¬ 
frared divergence of the ‘bare’ Coulomb interaction. Fol¬ 
lowing Ref.l^, we can write down the self-energy as a sum 
of line and residue term: 


E«(a;,k) = Er(a;,k) + Sr(cc,k), 

Ei‘“(a;,k) = -i^ f ^GAu^ + ^W 


(B3a) 




E(zfi,q)_ 


Fss'iKk'), 


(w,k) = - ©(-Cs'k')] 


qs' 

X 




A^s'k' - w,q)_ 


(B3b) 


F,,,(k,k'), (B3c) 


where E*/"® is obtained by carrying out the analytic con¬ 
tinuation iu>n —> w-|-*0 in Eq. (Bl) before performing the 
Matsubara sum and the residue term gives the correction 
to the self-energy due to the noncommutativity of these 
two operations. Using dimensionless energy (E = Ti/Ep, 
V = uj/Ep, Q = n/Ep) and dimensionless momenta 
(x = k/kp, x' = k'/kp, y = q/kp), we can write the 
self-energy as 


si™(ux) = 
sr(ux) = 




—a 




dy 


dy 


d^d 

(27r)^ 

d^d 


OdFss/(x,x') J 


dCl 


27r e(jD,y) p + — s'|x'| -|- 1 


(B4a) 


' e(»'|x'|-l-‘',y) 


(B4b) 
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where a = /(kvf) and 51^ is the solid angle element of dimension d. 

The chemical potential is determined by the electron energy at k = kp, i-e., 

^ = -^ + E.(0,l). (B5) 

Jl/p H/p 

For the self-energy at = 0 and x = 1, the residue term is zero and hence Ss(0,1) = Es"®(0,1). So we are left with 
evaluating only 1). The calculation of the self-energy is regulated by using momentum and frequency cutoff 

£ = kc/kp-. 


Ss(0,l) 



dCt 1 r d£ld -Fss/(l,x') 

27r J {2pY — s'|x'|-|-1 


(B6) 


Since the ultraviolet divergence stems from the integration over high-momentum region, we introduce the following 
parametrization 


y = Psmx, (B7a) 

Cl = Pcosx, (B7b) 

where 0 < % < tt, and evaluate the divergence part of the self-energy by focusing on the large-P region. In this region, 
we have 


f d£ld „ F'ss'(l)X') 

■“7 {2-kY — s'jx'l-I-1 


{ -^cosx + 
-^cosx + 


7rP2 L 
1 

152 


cos2x- -(2 + cos2x) 


O 


[cos 2x - s cos^ ^ 



(3D), 

(2D). 


(B8) 


Since the term of order 0{1/P) is odd in cosx, it will 
be zero after doing the integral over x- It remains then 
to evaluate the term C>(l/P^). Since we are interested 
only in the divergence part of the self-energy, in the self¬ 
energy calculation we use the intrinsic dielectric function 
which is given by 


e(*0,y) 


where 


1 + ^ln 


1 + 


GtT y2 + 
gpa y 


y2 + D2 


P 

1 + Aln- 
= 1 -I- A sin X 


(3D), 

(2D), 

(B9) 


A = 



(3D), 

(2D). 


(BIO) 


The dielectric functio n for 3D and 2D Dirac materials are 
given in Refs.mEsEZl and Refs.f^SITSl, respectively. Using 
Eqs. (B 81 and (B9), we can then calculate the self-energy 
as 


7-('-£-) 

[/o(A) + s/i(A)]ln£ 


(3D) 

-I-finite, (Bll) 

(2D) 


where 


/o(A) 


and 

/i(A) 



cos 2x 
1 + A sin X 
TT 2 — A^ 

A ~ 


arccos A 


/>7r 9 

a / cos X 

27r Jq ^1 + Asinx 


a 

ttA 



yr^A 2 

A 


arccos A 


(B12a) 


(B12b) 


The above derivation is given in Ref.l^. The results ob¬ 
tained above are consistent with those given in Refs.^^^^f^ 
and Ref.l^ for 2D and 3D Dirac materials, respectively. 

Considering extrinsic Dirac materials with Ep > 0 
(s = -l-l), we then obtain the zero-temperature RPA in¬ 
compressibility (up to the divergent part of the interact¬ 
ing term) as 


Ko 
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37r kp J 
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(B13) 
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where / 2 (A) = ^ (l - ^ + ^rccos a) . For 

the purpose of comparison, below we give the zero- 
temperature HF incompressibility (up to the leading- 
order divergent term in the interacting part) for 3D and 
2D Dirac materials: 


Ko 

K 


2a, 

1+“in 
4 


kr 

'kk 

(3D), 

kc 

kp 

(2D). 


(B14) 


FIG. 8. (Color online) (Top) Zero-temperature incompress¬ 
ibility of 3D Dirac materials as a function of wave vector cutoff 
kc/kp for different values of effective fine-structure constant: 
(a) a = 2.2 and (b) 0.01. (Bottom) Zero-temperature incom¬ 
pressibility of 2D Dirac materials as a function of wave vector 
cutoff kc/kp for (c) a = 2.2 and (d) 0.01. The above plots 
are calculated using the degeneracy factor g — 4. 

— HFA (kjkf = 10^) — HFA (k,/kf = W*) — HFA (kjkf = W-) 

— RPA (kjkf = 10'’) - RPA (kjk,, = lO"*) — RPA (kjkf = 10^) 




FIG. 9. (Color online) Zero-temperature incompressibility as 
a function of a for several different values of k^/kp of (a) 
3D Dirac and (b) 2D Dirac materials. Inset: Plot of zero- 
temperature incompressibility for small values of a. The 
above plots are calculated using the degeneracy factor g = 4. 


So within the RPA, the incompressibility Kq/K for 
a 3D Dirac material increases slowly with kc/kp with a 
log — log dependence while Kq/K for a 2D Dirac material 
increases as h\{kc/kp). Within the HFA, the incompress¬ 
ibilities Kq/K for 2D and 3D Dirac materials increase 
as \n{kc/kp). The dependence of the zero-temperature 
HF and RPA incompressibilities on the momentum cut¬ 
off kc/kp is shown in Fig.[^ Figure]^ shows the plots of 
zero-temperature incompressibilities as functions of the 
effective fine-structure constant a. As shown in the fig¬ 
ure, the HF incompressibility increases linearly with a 
while the RPA compressibility increases with a rate that 
becomes smaller as a increases where the exact depen¬ 
dence given in Eq. (B13). Note that for a ^ 1, the 
leading-order divergent piece of the HF and RPA com¬ 
pressibilities are the same. This can be seen from the 
slope of the plots depicted in the right panel of Fig. 
and insets of Fig. The fact that the leading-order in¬ 
teracting results in the coupling constant (for a ^ 1) 
of the compressibility agree for RPA and HFA indicates 
that we have no way of ascertaining which one is a better 
approximation within the single-loop theory, and exper¬ 
imental results are necessary in order to make progress. 
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